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An Analysis of the Weak Finite Element Method for Convection-Diffusion 

Equations 

Tie Zhang and Yanli Chen 

Department of Mathematics and the State Key Laboratory of Synthetical Automation for 
Process Industries, Northeastern University, Shenyang 110004, China 

We study the weak finite element method solving convection-diffusion equations. A weak 
finite element scheme is presented based on a spacial variational form. We established a 
weak embedding inequality that is very useful in the weak finite element analysis. The 
optimal order error estimates are derived in the discrete df^-norm, the L 2 -norm and the 
Loo-norm, respectively. In particular, the df^-superconvergence of order A; -|- 2 is given 
under certain condition. Finally, numerical examples are provided to illustrate our the¬ 
oretical analysis. 
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1 Introduction 

Recently, the weak Galerkin finite element method attracts much attention in the field of 
numerical partial differential equations [D El El H Em El 0 [9]. This method is presented 
originally by Wang and Ye for solving general elliptic problems in multi-dimensional do¬ 
main [T]. Since then, some modified weak Galerkin methods have also been studied, for 
example, see [iniEIlEllE3]. In general, a weak finite element method can be considered 
as an extension of the standard finite element method where classical derivatives are 
replaced in the variational equation by the weakly defined derivatives on discontinuous 
functions. The main feature of this method is that it allows the use of totally discontinu¬ 
ous finite element function and the trace of finite element function on element edge may 
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be independent with its value in the interior of element. This feature make this method 
possess all advantages of the usual discontinuous Galerkin (DG) finite element method 
[made], and it has higher flexibility than the DG method. The readers are referred 
to articles [HEllIl] for more detailed explanation of this method and its relation with 
other finite element methods. 

In this paper, we study the weak finite element method for convection-diffusion equa¬ 
tions: 


— div{A'Vu) + b ■ Vu + cu = f, in Q, (1.1) 

where D C coefficient matrix A = {aij)dxd and b is a vector function. 

We first introduce the weak gradient and discrete weak gradient following the way in 
[T]. Then, we consider how to discretize problem (jl.lh by using weak finite elements. In 
order to make the weak finite element equation have a positive property, we present a 
spacial weak form for problem (II.Ih . This weak form is different from the conventional 
one and is very suitable for the weak finite element discretization. We establish a discrete 
embedding inequality on the weak finite element space which provides a useful tool for 
the weak finite element analysis. In analogy to the usual finite element research, we 
derive the optimal order error estimates in the discrete fl^-norm, the L 2 -norm and the 
Loo-norm, respectively. In particular, for the pure elliptic problems in divergence form 
(6 = 0, c = 0), we obtain an 0(/i^"*“^)-order superconvergence estimate for the gradient 
approximation of the weak finite element solution, when the {k, A:-(-l)-order hnite element 
polynomial pair (interior and edge of element) is used. Both our theoretical analysis and 
numerical experiment show that this weak finite element method is a high accuracy and 
efficiency numerical method. 

This paper is organized as follows. In Section 2, we establish the weak finite element 
method for problem (|l.ll) . In Section 3, some approximation functions are given and 
the stability of the weak finite element solution is analyzed. Section 4 is devoted to the 
optimal error estimate and super convergence estimate in various norms. In Section 5, 
we discuss how to solve the weak finite element discrete system of equations and then 
provide some numerical examples to illustrate our theoretical analysis. 

Throughout this paper, for a real s, we adopt the notations IT*’P(D) to indicate the 
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usual Sobolev spaces on domain D equipped with the norm || • ||s,p,D and semi-norm 

1 • and if p = 2, we set = H^{D), || ■ ||s,p,D = || • IUd- When D = 12, we 

omit the index D. The notations (•, •) and || • || denote the inner product and norm in the 
space L2(12), respectively. We will use letter C to represent a generic positive constant, 
independent of the mesh size h. 

2 Problem and its weak finite element approximation 

Consider the convection-diffusion equations: 

{ —div(AVu) -|- b ■ Vu + cu = f, in 
u = g, on 912, 

where 12 C {d = 2,3) is a polygonal or polyhedral domain with boundary 912, coeffi¬ 
cient matrix A = {aij)(ixd is uniformly positive definite in 12, i.e., there exists a positive 
constant oq such that 

Ve G X € 12. (2.2) 

We assume that aij{x) G [fT^’°°(12)]'^^'^, h G [fT^’°°(12)]'^ and c(x) G Loo(ri). As usual, 
we further assume that 

c(x) — ^div6(x) >0, X G 12. (2.3) 

Let T/j = U{-^} bs S’ regular triangulation of domain 12 so that 12 = 
where the mesh size h = maxhx, hx is the diameter of element K. 

In order to define the weak finite element approximation to problem (12.ip . we first 
need to introduce the concepts of weak derivative and discrete weak derivative, which 
are originally presented in [T]. 

Let element A G T/i, denote the interior and boundary of K by and dK, re¬ 
spectively. A weak function on element K refers to a function v = with 

€ L 2 {K^), = v\qk G L 2 {dK). Note that for a weak function v = 

may not be necessarily the trace of on dK. 

Definition 2.1. The weak derivative d^.v of a weak function v with respect to vari¬ 
able Xi is defined as a linear functional in the dual space whose action on each 
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q € H^{K) is given by 

f d'^i'vqdx = — f v^dxiQdKdx + f v’^qcosOids, Vg G H^{K), 1 < i < d. (2.4) 
Jk ‘ Jk JdK 

where n = (cos 9i,..., cos 0^)^ is the outward unit normal vector on dK. 

Obviously, as a bounded linear functional on d^.v is well defined. Moreover, 

for V G H^{K), if we consider n as a weak function with components = v\xo, = v\qk, 

then by the Green formula, we have for q G H^{K) that 

[ dxiVqdx = - [ vdx,qdx+ [ vq cos 9ids = - [ v^dx,qdx+ [ v\ cos 9ids, (2.5) 
Jk Jk JdK Jk Jok 

which implies that d'^.v = dxiV is the usual derivative of function v G H^{K). 

According to Definition 2.1, the weak gradient of a weak function v should be 

such that VwV = {dx^v,... G satisfies 

f VwV-qdx = — f v^divqdx+ f v^q ■ nds,y q & (2.6) 

Jk Jk JdK 

Below we introduce the discrete weak gradient which is actually used in our analysis. 
For any non-negative integer I > 0, let Pi{D) be the space composed of all polynomials 
on set D with degree no more than 1. Introduce the discrete weak function space on K 


W{K-, k,r) = {v = G Pk{K^), G Pr{dK) }. 


(2.7) 


Definition 2.2. For v G W{K,k,r), the discrete weak derivative d^.^v G Pr{K) is 
defined as the unique solution of equation: 

I djf. ^vqdx = — I v^qdx -|- [ v^qcos9ids, Vq G Pr{K), 1 <i <d. (2.8) 

Jk ’’’ Jk JdK 

According to Definition 2.2, for weak function v G W{K,k,r), its discrete weak 
gradient = {dx^ ^'’^,... G [Pr{KY is the unique solution of equation: 

/ '^w,rV-qdx = — / v^dwqdx + / • nds, Vq G [Pr(A)]'^. (2.9) 

Jk ’ Jk JdK 

Remark 2.1. IFe here first define the (discrete) weak derivative, and then the (dis¬ 
crete) weak gradient follows naturally. This method of defining (discrete) weak gradient 
is slightly different from that in w, in which the (discrete) weak gradient is defined solely. 
From (12.6p and (|2.9|1 . we have 
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This shows that V-w,rV is a discrete approximation of in [Pr{K)]'^. In particular, if 
V € H^{K), we have from (I2.5jl and (12.9p that 


[ {Vv - ■ qdx = 0, Vq G [Pt{K)Y■ 


JK 

That is, Vw,rV is the L 2 projection of Vn in [Pr{K)]'^ if u G H^{K). 

We have introduced the weak derivative (gradient) and discrete weak derivative (gra¬ 
dient), but only the discrete weak gradient given in (j2.9p will be used throughout this 
paper. The others also should be useful in the study of numerical partial differential 
equations. 

A important property of can be stated as follows, see [H Lemma 5.1]. 

Lemma 2.1. Let v = G W{K,k,r) be a weak function and r > k. Then, 

'^w,rV = 0 on K if and only if v = constant, that is, = constant on K. 

Now, we construct the weak finite element space. Denote by = |J{ e G dK \ dkl : 
K G Th\ the union of all boundary faces or edges {d = 2) of elements in that are 
not contained in dLl. Let Ki and K2 be two adjacent elements with the common face 
e = dKi Pi dK 2 , and ni and 112 are the outward unit normal vectors on dKi and dK 2 , 
respectively. For weak function defined on Th, set v\ki = v\k 2 = {^ 2 )^ 2 }- We 

define the jump of weak function u on e by 

[v]e = v\ni + vln 2 = {v\ - V 2 )ni, e G 


Then, weak function v is single value on e if and only if [v\e = 0. The weak finite element 
space is now defined by 

Sh = {v. v\k G W{K, k, r), K G n, [uje = 0, e G 

S’)] = {u : V £ Sh, ujan = 0}. 


In order to define the weak finite element approximation to problem (|2.1I) . we need 
to derive a spacial weak form for problem (12.ip . From the differential formula 

b • Vtt = -b ■ Vu + -div(bM) — -divbu, (2T0) 

and the Green formula, we see that a weak form for problem (12.ip is to find u G FI^(D) 
u\dn = 9 such that 

(AVu, Vu) -b i(& • Vu,u) - ^{u,b- Vu) -t- {cbU,v) = (/,u), Vu G (2.11) 
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where Cb = c — ^div6 > 0. Denote the discrete L 2 inner product and norm by 

{u,v),= ^iu,v)K= Y: [ --dx, \\u\\l = iu,u),. 

KdTn K&Th 

Motivated by weak form (12.lip , we define the weak finite element approximation of 
problem (j2.ip by finding Uh € S^, Uh\dfi = Qh such that 

ah{uh,v) = (/,u°), Vu G S’”, (2.12) 

where gh is a proper approximation of function g and the bilinear form 

u) = (AVuijrU, ^'w,rV)h ^ w,r'lii V )}i 2 ^ ipbU , V ). (2.13) 

Bilinear form ah{u,v) is not based on the conventional one: 

a{u, v) = (AVu, Vu) + {b ■ Vu, v) + {cu, v). 

The advantage of our bilinear form is that it always is positive definite on the weak 
function space x Sh, and the conventional one is not, since the integration by parts 
does not hold on weak function space Sh or S'”. 

Theorem 2.1. Let r > k. Then, the solution of weak finite element equation 112 . 12 \) 
uniquely exists. 

Proof. Since equation (12.121) is essentially a linear system of equations, we only 
need to prove the uniqueness. Let f = gh = 0, we need to prove Uh = 0. Taking 
u = G 5” in (I2.12p . we obtain 

UQIIII/j ^ ah{Uh,Uh') — 0. 

This implies that ^w,rUh = 0 on T/j. Thus, from Lemma 2.1, we know that Uh is a 
piecewise constant on Th. Since [uh]e = 0 and Uh\dn = 0, so we have Uh = 0. ■ 

3 Projection and approximation 

In this section, we give some projections and approximation properties which will be 
used in next section. 

In order to balance the approximation accuracy between spaces Sh and Pr{K) used 
to compute S/w,rUh, from now on, we always set the index r = A: + 1 in (I2.7l) - (l2.8p . The 
other choice of weak finite element space can be found in mm- 
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(3.1) 


(3.2) 


(3.3) 


For I > 0, let is the local L 2 projection operator, restricted on each element K, 
Pj^ : u L 2 {K) —>■ G Pi{K) such that 

{u- Plu,q)K = 0, Vg G Pi{K), K G T/j. 

By the Bramble-Hilbert lemma, it is easy to prove that (see m) 

||u — Pf^u\\Q^x ^ ) 0 < S < / + 1. 

We now define a projection operator Qh '■ u G 77^(0) QhU € Sh such that 

Qhu\K = {Qlu,Qiu] = {P^u,Pl+^v!>], K G n, 

where Pq^^ is the L 2 projection operator in space Pk+i{dK). 

Lemma S.l.E Let u G s > 0. Then, QhU has the following approximation 

properties 

11'^ “ ^ 0<s<A: + l, K €z Th, (3-4) 

\\"^w,rQhU - Vu\\o,K < Ch''j^\\u\\i+s,K, 0 < s < k + 2, K eTh. (3.5) 

Proof. Since Q^u = P^u, then estimate ()3.4p follows from (j3.2p . Furthermore, 
from ()2.9p and the definition of QhU, we have 

/ ^w,rQhU ■ qdx = — Q^udivqdx + / Q\uq ■ nds 

Jk ’ Jk JdK 

= — udivqdx + / uq ■ nds = / Vttqdx, Vq G [Pr(A)]^. 

Jk JdK Jk 

This implies Vw,rQhU = Pff^u and estimate (13.5p holds, noting that r = k + 1. ■ 

For the error analysis, we still need to introduce a special projection function m- 
For simplifying, we only consider the case of two-dimensional domain {d = 2). 

Let Ci and Aj (1 < i < 3) are the edge and barycenter coordinate of K, respectively. 
For function cj), curie/) = {dx2<f, —dxi 4 >)'^- Let Space P^_|_2(A) = {p G Pk+2{K) '■ p\dK = 
0} = AiA 2 A 3 Pfc_i(A) and H{div,Q) = {u G [L2(A1)]^ : divu G L 2 {Ll)}. 

Define the projection operator tt/i : 77(div;D) ^ 77(div;D), restricted on A G T^, 
TT/iit G [Pk+i{K)]^ satisfies 


{u - ■KhU,Vq)K = 0, Vg G Pfc(A), 

{u- TThu) ■ nqds = 0, Vg G Pfc+i(ej), i = 1,2,3, 

J ei 

{u - TThU, curl q)K = 0, Vg G P°+ 2 (A). 


(3.6) 

(3.7) 

(3.8) 
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Some properties of projection Tr^u had been discussed in we here give a more detailed 
analysis for our argument requirement. 

Lemma 3.2. For u € H{div;Q), the projection -KhU uniquely exists and satisfies 

{div{u - 7rhu),q)K =0, Vg € Pk{K), K € Th- (3.9) 

Furthermore, if u G s > 0, then 

\\'^hu\\o,K < C'||ii||i,x, K G Th, ( 3 . 10 ) 

\\u — TThuWo^K < Ch%;\\u\\s,K, I < S < k + 2, K £Th. ( 3 - 11 ) 

Proof. We first prove the unique existence of Tr^rt. Since the number of dimensions 
(noting that (13.611 is trivial for q = constant)'. 

dim{Pk{K)) - 1 + 3dim(Pfc+i(ej)) + dim(P^+ 2 (A)) = 2dim(Pfc+i(A)), 

so the linear system of equations (I3.6p ~ (l3.8p is consistent. Thus, we only need to prove 
the uniqueness. Assume that rt = 0 in ()3.6P '^ (j3.8p . we need to prove iThU = 0. From 
dMP-dlTl), we have iThU • n = 0 on dK and 

(div7r/jU,g)K =-(vTfeU, Vg)x + / tthU ■ nqds = 0, y q e PkiK). 

JdK 

This implies divvr/iU = 0 on A. So there exists a function cj) G Pk+ 2 {K) so that curl(/> = 
TTfiU (see m)' Since the tangential derivative drfi = TThU • n = 0 on dK, so = fio = 
constant on dK. Let p = — (j)Q. Then, p G P^_|_ 2 (A) and curlp = iThU. Taking 

q = p in (13.8p . we obtain = 0 so that tt/jU = 0. Next, we prove (13.911 - 

(j3.1ip . Equation (13.9p comes directly from the Green formula and (|3.6p - (l3.7p . From the 
solution representation of linear system of equations (I3.6p ^ ()3.8p . it is easy to see that 
on the reference element K, 

P 11 ^ 110 , 911 ) — ll^^llo,.;!)’ (3-12) 


where we have used the trace inequality. Then, (I3.10p follows from (|3.12l) and a scale 
argument between K and K. From (I3.12p . we also obtain 
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Hence, estimate (I3.1ip can be derived by using the Bramble-Hilbert lemma. ■ 

The following discrete embedding inequality is an analogy of the Poincare inequality 
in Hq(H). 

Lemma 3.3. Let Q be a polygonal or polyhedral domain. Then, for v € S^, there is a 
positive constant Cq independent of h such that 

lb°|| < Co||V^,,u|U, Vu € Si (3.13) 

Proof. For v e S“. we first make a smooth domain D H ( if H is convex, we 
may set = H) and extend to domain Q' by setting = 0. Then, there exists 

a function w G Hq{Q') f) such that 

—Aw = in Q', ||rc|| 2 ,n' < C'||u'^||. 

Now we set w = —Vw, then w € satisfies 

divm = in Ll, ||rr||i < ||ru|| 2 ,n' < C'||u°||. 


Hence, we have from (13.9p . (|3.10p and (|2.9p that 

||u°|p = (divm,u°) = (divTr/ite, u°) 


'^w,rV ■ TThWdx + / V^TTfiW ■ nds) 

rri. Jk ’ JdK 




^ / Vu,,rV -T^hWdx < \\S/w^rV\\h\\'n'h'w\\h 

Ken 


IK 


< C'||V^,,.u||h||m||i < CllV^^r-ulIhllu^ll, 


where we have used the fact that [u]e = 0 and 



(3.14) 


The proof is completed. ■ 

A direct application of Lemma 3.3 is the stability estimate of weak finite element 
solution Uh- 

Lemma 3.4. Let Uh G Sh be the solution of problem ph = Q\g and u G H^{LL) 

the solution of problem 112.1\) with f G L 2 {Ll) and g G H^{dLL). Then we have 


ll4ll + llv.,.«/^IU<c(||/|| + |bl|l,9^)• 
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Proof. Let = Uh — QhU- Then, from (|2.12p . we see that eh G satisfies 
ah{eh,v) = {f,v^) - ah{QhU,v), v G 8°- 
Taking v = eh and noting that Q^u = Pj^u and '\/w,rQhU = Ph'^'^j we have 

||Vi„,reh||h — 11/1111^/111 + (l^loo + \b\oo + |c6|oo)(||Vu|| + ||t(||) (|| Vi„^re/i ||/i + ||e°||). 

Using embedding inequality (|3.13p and the a priori estimate for elliptic problem (|2.ip . 
the proof is completed. ■ 

4 Error analysis 

In this section, we do the error analysis for the weak finite element method (I2.12p . We 
will see that the weak finite element method possesses the same (or better) convergence 
order as that of the conventional finite element method. 

In following error analysis, we always assume that the data A, b and c in problem 
m is smooth enough for our argument. 

Lemma 4.1. Let u G be the solution of problem 112.1\) . Then we have 

{-j:h{AVu),Vw,rV)h + • Vn, - ]^{TTh{hu), Vy,^rV)h + (Cfett, u°) = (/, u°), V G Sl- 

Proof. Let w G From Lemma 3.2 and (I3.14p . we have 

- (divr(;,u°)/i = -{diYTrhW,v^)h = {TrhW,\/^^rv)h, v G S'}]. (4.1) 

Next, from equations (|2.ip and (|2.10l) . we obtain 

-(div(AVu), u°)fe + ^ib- Vu, v^)h - ^(div(bn), u°)h + {cbU, u°) = (/, u°), v G S^. 

Together with (14.ip in which setting w = A'Vu and w = bu, respectively, we arrive at 
the conclusion. ■ 

We first give an abstract error estimate for Uh — Qh'U' in the discrete Ll^-norm. 
Theorem 4.1. Letu anduh be the solutions of problems f2.1\} and f2.12\} . respectively, 
u G and ph = Qh9- Then, we have 

O‘0\\^w,r‘Uh ^w,rQh^\\h — \\AS/w,rQh'^ ^/i(4V'u) ||/i + C*o||fj • {^w,rQh^ ^^)||/i 

+ \\bQlu - Trh{bu)\\h + Co\\cb{QhU - u)\\h. 


(4.2) 
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Proof. From Lemma 4.1, we see that QhU satisfies the equation 
ah{QhU,v) = {f,v^) + {AS/w^rQhU- 7rhiAS/u),Vw,rv)h 

^ ^ )h 2 w h 

+{cb{Qlu - u),v°)h, Vv e 8°. (4.3) 

Combining this with equation (|2.12p . we obtain the error equation 

^ )/i- 2 ^w,r^)h 

+(cb(QhU - u),v°)h, Vv € 8°. (4.4) 

Taking v = QhU — Uh in (14.41) and nsing embedding inequality (I3.13p . we arrive at the 
conclusion of Theorem 4.1. ■ 

By means of Theorem 4.1, we can derive the following error estimates. 

Theorem 4.2. Letu anduh be the solutions of problems 112.1\} and 112.1^) . respectively, 
u G s > 0, and gh = Q\g- Then we have the optimal order error estimates 

\\uh - "“11 + \\'^w,rUh - < Ch^^^\\u\\ 2 +s, 0 < s < k. (4.5) 

In particular, for the pure elliptic problem in divergence form (b = 0,c = 0) and u G 
s > 0, we have the superconvergence estimate 

\\V^,rUh-Vu\\h<Ch^^^\\u\\ 3 +s) 0 ^ s ^ k. (4-6) 

Proof. Using the approximation properties (j3.4p ~ f|3.5p and (|3.1ip . we obtain 

\\^^w,rQh'^ II 

< \A\oo\\'^w,rQhU - Vn||/i + ||AVtt - Trh{AVu)\\h < Ch^~^^\\u\\ 2 +s, 0 < s < k + 1, 
\\b-{S/^,rQhU-S/u)\\h<Ch^+^\\u\\ 2+s) 0 < s < A: + 1, 

\\bQlu - 7rh{bu)\\h + ||cfe(Q°u - u)\\h < C'/i^+^||rt||i+s, 0 < s < A;. 

Substituting these estimates into (14.211 , we obtain 


^w,rQh‘^\\h — Ch || n || 2 + 5 , 81 S 81 k . 


(4.7) 
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Hence, using the triangle inequality 


^u},rUh ^ ^w,rQh'^\\ “ 1 “ || 


and approximation property (13.5p . estimate (14.5p is derived for the discrete ff^-norm. 
Since 


4 - < ||4 - Qlu\\ + WQlu - u 


then the L 2 -error estimate follows from the discrete embedding inequality (13.131) . esti¬ 
mates (I47fll and (lOI) . 

Furthermore, if 5 = 0, c = 0, from Theorem 4.1, we have 

flo||^ui,r^h 'w,rQh'^^h ^ ||/i' 

Then, the superconvergence estimate (|4.6I) can be derived by using approximation prop¬ 
erties ()3.5p and (13.lip and the triangle inequality. ■ 

Theorem 4.2 shows that the weak finite element method is a high accuracy numerical 
method, in particular, in the gradient approximation. 

From Theory 4.2, we see that, for A:-order finite element, weak finite element method 
usually has higher accuracy than other finite element methods in the gradient approxi¬ 
mation. The reason is that the discrete weak gradient is computed by using (A:-|-1)-order 
polynomial. This action will add some computation expense, but all additional compu¬ 
tations are implemented locally, in the element level, see (|2.9I) and Section 5. 

Below we give a superclose estimate for error Q^u — u^. To this end, we assume prob¬ 
lem (j2.ip has the Ff^-regularity and consider the auxiliary problem: w € FFg (H) P| 
satisfies 


— div(AVr(;) — b ■ Vw + cw = Q^u — Uh, in ||rc ||2 < C\\Q^u — u^\\. (4.8) 


From the argument of Lemma 4.2, we know that w satisfies equation: 




(4.9) 


Theorem 4.3. Letu anduh be the solutions of problems 112.1\} and i2.12\) . respectively, 
u G >0,gh = Q\g- Then we have the following superclose estimate 


WQ^u - u°|| < C'/i^+^||u|| 2 +s, 0 < s < k 


(4.10) 
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and the optimal L 2 -error estimate 

\\u^ — ^i|| < Ch^~^^\\u\\k+i, k > 1. (4-11) 

Proof. Taking v = QhU — Uh in (14.91) . we have 

WQhU - ulf = (V^^riQhU - Uh),7rh{AVw))h 
1 1 

--{b ■ Vw, Qlu - u\)h + -{T^h{bw),VwAQhU - Uh))h + {cbiQh'^ “ ul),w) 

= C^w,r{QhU - Uh),7rh{AVw) - AVy,^rQhw)h {Vw - Vw,rQhw),Qlu - ul)h 

+^{TTh{bw) - bQlw, Vw,r{QhU - Uh))h + {Cb{w - QhW), QhU - Uh) 

+ah{QhU - Uh,Qhw) 

^ (i7/l||'in||2(||V w,r{QhU Uh) || /i T || II) T {QhU Uh iQhw) 

< Ch2+"||n|| 2 +s||w||2 + ah{QhU - Uh, Qhw), 0 < s < A:, (4.12) 

where we have used embedding inequality (13.131) . estimate (14.71) and the approxima¬ 
tion properties of -KhW, V^.rQhW and Q^w. Below we only need to estimate ah{QhU — 
Uh,Qhw)- Using error equation (14.4p . we have 

^h{QhU Uh,,QhU!) — {AV^^rQhU Trh{AVu'),Vm^rQhUj)h 

+ ^{b ■ {^w,rQhU - Vu), Qhw)h - ^{bQhU - vr/i(bn), Vy,^rQhw)h 
+ {cb{Qh,u - u), Qhw)h, Vn € 5° 

= El -|- E 2 + E 3 -|- E 4 . (4.13) 

Since Qh = and V^^rQh = PJ), then by using Green’s formula and Lemma 3.2, we 
have 


El = {AVw,rQhU - 7rh{AVu),V^u,rQhW -Vw)h + {AV^u,rQhU - TTh{AVu),Vw)h 
< C'/i^+'^||n|| 2 +<i||w '||2 + {AVyj^rQhU - AVu, Vw)h + {AVu - 7rh{AVu), Vw)h 
= Ch2+^||n||2+.|k||2 + {^w,rQhU - Vu,A^Vw - Pj^{A^Vw))h 

— (div(AVtt — ■Kh{AVu)),w — Phw)h < C'/i^'''^||n|| 2 +s||rc|| 2 , 0 < s < k. 
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Similarly, we obtain 

E 2 — w,rQhU V It) ^ w)fi 

+ ^{V^,^rQhU-'^u,bw - pjl{bw))h < Ch‘^~^’^\\u\\ 2 +s\\w\\i, 0 < s <k, 

E 3 = -^{bQlu-7^h{bu),V^,^rQhW-Vw)h-^{bQlu-Trh{bu),Vw)h 

< Ch^+^\\u\\2+s\\w\\2 - ^{bQlu- bu,Vw)h - ^{bu - Trh{bu),Vw)h 

= Ch^+^\\u\\2+s\\w\\2 - ^{QU - u, bVw - Pt{bVw))h 

< Ch‘^~^''\\u\\2+s\\w\\2, 0 < S < k, 

Ea = {cb{Qlu-u),Ql^w)h = {cb{Q\u-u),Q\w-w)h 

+{{QhU - u),CbW - Ph{cbw))h < C'/i^+*||n||i+s||t(;||i, 0 <s <k. 

Substituting estimates Ei ~ Ea into (j4.13p and combining (I4.12p . we arrive at estimate 
(I4.10p . noting that ||tc ||2 < C\\Q^u — u°||. Estimate (|4.1ip follows from (|4.10p and the 
triangle inequality. ■ 

The difference between estimates (14.5p and (j4.1ip is that for getting the 0{h^^^)- 
order error estimate in the L 2 -norm, the regularity requirement in (|4.11l) is optimal and 
lower than that in (14.5p . 

In order to derive the Loo-error estimate, we need to impose the quasi-uniform con¬ 
dition on partition T/j so that the finite element inverse inequality holds in Sh- 
Theorem 4.4. Assume that partition is quasi-uniform, and u and Uh are the solu¬ 
tion of problems 112.1\) and /i2.12\) . respectively, u G P| s > 0, gh = 

Q\^g. Then, we have 

Ik - ^^°||o,oo < + lkl|2-Hs)i 0 < S < /C. (4.14) 

Proof. From Theorem 4.3 and the finite element inverse inequality, we have that 
WQhU - ^^hllo,oo < Ch~i WQbU - n° II < Ikk-Hs- 

Hence, by using the approximation property of Q^u = Pj^u, we obtain 

Ik ^/illo,oo ^ Ik Qh'^\\o,oo T IIQh^ '^^/illOjOO 

< C/i^+^lklll-Hs.oo + Ch‘^'^^~^\u\\2+s- 
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The proof is completed. ■ 

For two-dimensional problem (d = 2), Theorem 4.4 gives the optimal order error 
estimate in the Loo-norm. 

5 Numerical experiment 

In this section, we discuss how to solve the weak finite element equation ()2.12p and give 
some numerical examples to illustrate our theoretical analysis. 

5.1 Weak finite element linear system of equations 

In order to form the discrete linear system of equations from weak finite element equation 
(j2.12[l . we first introduce the basis functions of space 5°. Let A is a element and e is an 
edge of K. Further let {ipj^K{x), j = 1,, Nj^} be the basis functions of space Pk{K), 
{(pj^e{x), j = 1,..., Nf>} the basis functions of space Pk+i{e) and {ipjfi{x), j = 1,..., Mh} 
be the basis functions of polynomial set {p € Pk+i{dK), K ^ : p\e = pj^eix), e € 

Set the weak basis functions '4’j,K = {v’j,K,0} and = {0,<Pj,b}- Then, we have 
Si = spanii/ji^K, ■ ■ ■ G Th}. By definition ([IJ]) of discrete 

weak gradient, we see that the support set of S/w,r'4’j,Kix) is in K and the support set of 
^w,r'4’j,b{x) is in = {K : dK{^supp{il^j^b) 7^ 0}- Thus, weak finite element equation 
(I2.12[) is equivalent to the following linear system of equations: Uh € Sh, Uh\dn = gn, 
such that 

aK{uh,v) = V G K G T/*, (5.1) 

aRbi'^h.v) = 0, u G i = 1, - • • (5.2) 

where aD{u,v) is the restriction of ah{u,v) on set D, i.e., all integrals in ah{u,v) are 
restricted on D. Equations (j5.ip - (j5.2l) form a linear system composed of A x Nk -I- Mb 
equations with N x + Mb unknowns, where N is the total number of elements in T^. 

To solve equations (I5.ip ~ (l5.2l) . we need to design a solver to compute the discrete 
weak gradient Vw,rV or Vw,rUh- According to (12.91) . for given v = G W{K,k,r) 

{r = k + I), Vyj^rV G [Pk+i{K)]'^ can be computed by the following formula 

MKV^^r = AkV^ + ^ 

eedK 


(5.3) 
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where 14 ,,r, and are the vectors associated with functions Vw,rV G [Pk+i{K)Y, 
G Pk{K) and = v^\e G Pfc+i(e), respectively. The matrixes in (15.3p are as follows 

Mk = {'mij)dNrXdNr! = {0'ij)dNrXNkJ Be = {bij)dNrXNe^ 

'kklij j,Ki,K)Kt ^ij b^j ' ^')e- 

Now, linear system of equations (I5.1I) ~ (I5.2I) can be solved in the following way. We first 
use formula (j5.3p to derive the linear relation SIw,rUh{K) = L{u\^{K),u\{dK)). Then, 
by substituting Vw,rUh{K) into equations (j5.iP " (j5.2|] . we can obtain a linear system of 
equations that only concerns unknowns {rt^(A), n^(9A)}, A G T^. This linear system 
can be solved by using a proper linear solver, in which Vy^^rv{K) is computed by formula 

dSS]). 

5.2 Numerical example 

Let us consider problem (12.ip with the following data: 

A = (1 + X 1 X 2 )/, b = (1, 2)"^, c = sin(xiX 2 ), u = sin(xi7r) sin(x 2 '/r), (5.4) 

and the source term / = —div(AVn) + b ■ S/u + cu, domain = (0,1)^. 

In the numerical experiments, we first partition into a regular triangle mesh 

with mesh size h = 1/N. Then, the refined mesh T/j /2 is obtained by connecting the 

midpoint of each edge of elements in Th by straight line. Thus, we obtain a mesh series 

= 0,1 ,... . We use the polynomial pair {Po{K), Pi{dK)) for space W{K,k,r) 

and set Pr{K) = Pi{K) in the weak finite element discretization (12.121) . We examine the 

computation error in the discrete A^-norm, the L 2 -norm and the Loo-norm, respectively. 

The numerical convergence rate is computed by using the formula r = ln(e/i/e^)/ln2, 

2 

where Ch is the computation error. Table I give the numerical results. We see that the 
convergence rates are consistent with or better than the theoretical prediction. Then, we 
further examine the super convergence of the weak finite element solution (see Theorem 
4.2). We take the data as in (|5.4p with 6 = 0 and c = 0. The desired A^-superconvergence 
is observed from the numerical results, see Table II. In particular, from Table I-II, we see 
that the numerical convergence rates in L 2 -norm and Loo-norm are also superconvergent, 
although we have no such conclusion in theory. 
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TABLE I History of convergence 



w,r'^h 

<1 

5 ^ 

\\u-ul 

II 

IK - 

'^lloo 

mesh h 

error 

rate 

error 

rate 

error 

rate 

1/4 

0.9329 

- 

0.870e-l 

- 

1.990e-l 

- 

1/8 

4.660e-l 

1.0013 

2.427e-2 

1.8417 

5.525e-2 

1.8489 

1/16 

2.325e-l 

1.0033 

6.246e-3 

1.9579 

1.413e-2 

1.9677 

1/32 

1.162e-l 

1.0011 

1.573e-3 

1.9893 

3.555e-3 

1.9905 

1/64 

5.806e-2 

1.0003 

3.940e-4 

1.9973 

8.906e-4 

1.9969 

1/128 

2.903e-2 

1.0001 

0.986e-4 

1.9993 

2.228e-4 

1.9993 


TABLE 11 

History of convergence with b = 

0 and c = 0 



II V-u;,rR/l “ 

< 

Ik - < 

II 

K- 

'^lloo 

mesh h 

error 

rate 

error 

rate 

error 

rate 

1/4 

1.875e-l 

- 

3.129e-2 

- 

0.928e-l 

- 

1/8 

4.896e-2 

1.9370 

8.538e-3 

1.8735 

2.553e-2 

1.8625 

1/16 

1.239e-2 

1.9821 

2.184e-3 

1.9673 

6.533e-3 

1.9664 

1/32 

3.109e-3 

1.9948 

5.490e-4 

1.9917 

1.642e-3 

1.9919 

1/64 

7.782e-4 

1.9984 

1.375e-4 

1.9979 

4.112e-4 

1.9979 

1/128 

1.946e-4 

1.9995 

3.437e-5 

1.9995 

1.028e-4 

1.9993 
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